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Magnetic and Nonmagnetic Phases in Doped AB 2 t-J Chains 


R. R. Montenegro-Filho and M. D. Coutinho-Filho 
Laboratorio de Fisica Teorica e Computacional, Departamento de Fisica, 
Universidade Federal de Pernambuco, 50670-901, Recife-PE, Brazil 

We discuss the rich phase diagram of doped AB 2 t-J chains using data from DMRG and exact 
diagonalization techniques. The J vs 5 (hole doping) phase diagram exhibits regions of itinerant 
ferrimagnetism, Incommensurate, RVB, and Nagaoka States, Phase Separation, and Luttinger Liq¬ 
uid (LL) Physics. Several features are highlighted, such as the modulated ferrimagnetic structure, 
the occurrence of Nagaoka spin polarons in the underdoped regime and small values of J = 4t 2 /[7, 
where t is the first-neighbor hopping amplitude and U is the on-site repulsive Coulomb interaction, 
incommensurate structures with nonzero magnetization, and the strong-coupling LL physics in the 
high-doped regime. We also verify that relevant findings are in agreement with the corresponding 
ones in the square and n-leg ladder lattices. In particular, we mention the instability of Nagaoka 
ferromagnetism against J and 8. 


I. INTRODUCTION 

The t-J version of the Hubbard Hamiltonian [l] is a 
key model for the understanding of strongly correlated 
electron systems. The model is defined through only 
two competing parameters: the hopping integral t , which 
measures the electron delocalization through the lattice, 
and the exchange coupling J = 4 t 2 /U, where U » t is 
the on-site Coulomb repulsion. In fact, several versions of 
the simplest Hubbard Hamiltonian, with a single orbital 
at each lattice and the on-site Coulomb repulsion, have 
been extensively used to model a variety of phenomena, 
such as: metal-insulator transition Q, quantum mag¬ 
netism Q and High-T c superconductivity [l- Moreover, 
exact solutions p] and rigorous results [Bj, ;6| have played 
a central role in this endeavor. 

We emphasize Lieb’s theorem jj], a generalization of 
the one by Lieb and Mattis [Ij for Heisenberg systems, 
which asserts that the ground state (GS) total spin of 
a bipartite lattice at half filling and U > 0 is given 
by Sgs = \Na ~ Nb |/2, where N ,4 ( Nb ) is the num¬ 
ber of sites on sublattice A ( B)\ indeed, Lieb’s theorem 
has greatly enhanced the investigation of new aspects of 
quantum magnetism In particular, we mention the 
occurrence of ferrimagnetic GS, in which case we select 
studies using Hubbard or t-J models LL 15|, including the 
Heisenberg strong-coupling limit [ 1 (ill 1 8} , on chains with 
AB 2 or ABC topological structures with Sgs = 1/2 per 
unit cell [3 1~H1 ilTil . [l7| , which implies ferromagnetic and 
antiferromagnetic long-range orders [To} . Further, the 
inclusion of competin g in teractions or geometrical and 
kinetic frustration |19i42I| . enlarge the classes of models, 
thereby allowing ground-states not obeying Lieb or Lieb 
and Mattis theorems. These studies have proved effective 
in describing magnetic and other physical properties of a 
variety of organic, organometallic, and inorganic quasi- 
one-dimensional compounds [l9|, [ 23 ] • 

Of particular physical interest are doped systems, al¬ 
though in this case rigorous results are much rare [f|. 
One exception is Nagaoka’s theorem 23|, which asserts 
that for J = 0 (U 00 ) the t-J model with one hole 


added to the undoped system (half-filled band) is a fully 
polarized ferromagnet , favored by the hole kinematics, if 
the lattice satisfies the so-called connectivity condition 
[24| . A long-standing problem about this issue is the sta¬ 
bility of the ferromagnetic state for finite hole densities 
and finite values of J. Numerical results have indicated 
[ 2 BL m that two-dimensional lattices display a fully po¬ 
larized GS for J = 0 and S < 0.2, where 8 = N^/N, with 
Nh (N) the total number of holes (sites); while, analyti¬ 
cal studies [13, [SI have suggested that this state is stable 
up to J t ~ S 2 . 

Further, an ubiquitous phenomenon in doped strongly 
correlated materials is the occurrence of inhomogeneous 
states, particularly spatial phase separation in nano- and 
mesoscopic scales fl3 | and incommensurate states [ 29 L 
M- In underdoped High-T c materials, dynamical and 
statical stripes in copper oxide planes has been the focus 
of intensive research [311 ] . Concerning two-dimensional t- 
J or Hubbard models, phase separation into hole-rich and 
no-hole regions was discussed in the large— and small— J 
limits [32|]. However the precise charge distribution in the 
ground state remains controversial. The use of distinct 
and refined numerical methods have pointed to striped 
[33l | or uniform phases [HJ ; recently, it was claimed that 
the origin of this issue relies on the strong competition 
between these phases [35}. For the linear t-J Hubbard 
chain the physics is more clear [36} , and phase separation 
takes place for J = 2.5 — 3.0, depending on the doping 
value, but it is absent in the small— J regime. 

In this work, we use Density Matrix Renormalization 
Group (DMRG) [37j technique and Lanczos exact diago¬ 
nalization (ED) to obtain the ground state phase diagram 
and the low-energy excitations properties of the doped t- 
J model on AB 2 chains [ 9 } for J = 0.0 — 0.4. We verify 
the occurrence of an itinerant modulated ferrimagnetic 
(FERRI) phase in the underdoped regime, regions of in¬ 
commensurate (IC) states and Nagaoka ferromagnetism 
(F), and two regions of phase separation (PS), in which 
IC and F states coexist with the resonating valence bond 
state (RVB), respectively. In addition, we find that the 
RVB state is the stable phase at 8 = 1/3, and identify 
a crossover region that ends at the onset of a Luttinger 
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FIG. 1. (Color online), (a) GS phase diagram for the AB 2 
t-J model (error bars account for the discrete values assumed 
by S in a finite-size system). The phases are illustrated in 
(b): modulated ferrimagnetism (FERRI), incommensurate 
(IC), Nagaoka ferromagnetism (F), short-range resonating va¬ 
lence bond (RVB) states, phase separation (PS), and Lut- 
tinger liquid (LL). The estimated transition lines 5ferri,j, 
Sps,j, and Jf,s are also pointed out. (c) Ground state total 
spin, Sgs, normalized by its value in the undoped regime: 
Sl = (N c / 2) — 0.5, as function of $ for the indicated values 
of J and N = 3 N c + 1 = 100. 


liquid (LL) phase at 6 = 2/3, above which the LL physics 
(38l | sets in. 


II. PHASE DIAGRAM 

The t-J model reads: 

= Pg(cIc j<7 + H.c.)P G (1) 

<i,j>,cr 

+ j ^2 (Si • S j - -TliTlj), 

<i,j> 

where Ci a annihilates electrons of spin a at site i, n, is the 
number operator at site i and Pq = IL(i - ni-friii) is the 
Gutzwiller projector operator that excludes states with 
doubly occupied sites. In our simulations, we set t = 1 
and have considered chains with N c (TV) unit cells (sites). 
In ED calculations closed boundary conditions are used 
with N c = 8 (N = 3 N c ), while in the DMRG simulations 
open boundary conditions are used and the system sizes 



FIG. 2. (Color online), (a) Effective linear chain (spacing 
a = 1) associated with N = 3 N c + 1 = 100 sites for J = 0.1 
used to illustrate the hole, ( rih,i ), and spin, (Si), profiles: (b) 
5 = 4/100 (FERRI phase) and (c) 5 = 18/100 (IC phase). 


ranged from N c = 33 ( N = 3TV C + 1 = 100) to N c = 121 
(N = 364). We retain from 243 to 364 states in the 
DMRG calculations, and the typical discarded weight is 
1 x lO" 7 . 

The ground state (GS) phase diagram, shown in Fig. |T| 
(a), displays the regions of the above-mentioned phases, 
illustrated in Fig. [2(b), including the estimated transi¬ 
tion lines and the crossover region. A special feature of 
the AB 2 chain is its symmetry [12|, [l3|, [l7| under the ex¬ 
change of the labels of the B sites in a given unit cell l 
[identified in the FERRI state, Fig. [2(b)]. This symme¬ 
try implies in a conserved parity pi = ±1 in each cell of 
the lattice. The phase diagram of a chain with N c unit 
cells is calculated by obtaining the lowest energy for all 
subspaces with x contiguous cells of parity —1 and the 
others N c — x cells with parity +1, with x = 0... N c , 
for fixed <5 and J. In the phase diagram shown in Fig. 
[2(a), p = Y^i=iPi = +1 f° r > 1/3, p / ± 1 in the PS 
region, and p = — 1 for S < Sps,j■ The magnetic configu¬ 
ration of a phase is identified by the total spin Sgs , local 
magnetization, magnetic structure factor, and spin cor¬ 
relation functions. In what follows, we shall characterize 
the phases shown in Fig. [2 a)- 


III. FERRIMAGNETISM AND TRANSITION 
TO IC STATES 


At S = 0 and J / 0, the insulating Lieb ferrimag- 
netic state with total spin quantum number Sgs = Sl = 
N c /2 — 0.5 = Sl is found for a chain with open bound- 
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FIG. 3. (Color online). Chain with N = 3N C + 1 = 100 and 
J = 0.3. (a) and (b): Magnetic structure factor S(q) for the 
indicated values of 8. Inset of (b): A q = g max — n, where q m ax 
is the value of q at which the local maximum of S(q), near 
q = 7r, is observed. 


ary conditions, N = 3 N c + 1 = 100, with an A site on 
each side. In order to evaluate the stability of this state 
against doping, we calculate Sgs as a function of 8 from 
the energy degeneracy in S z . As shown in Fig. H c ); 
as hole doping increases from 8 = 0 to a critical value 
8 = Sferri,j , the value of Sgs decreases linearly from 
Sl to 0 or a residual value, signaling a smooth transi¬ 
tion to the IC phase. However, for low enough J, Sgs of 
the IC phase increases linearly with 8 up to S = 8ps,j, 
the line at which PS occurs [see Fig. Oja)], or up to 
the boundary, Jf,s, of the Nagaoka F phase. This unex¬ 
pected behavior claims for an explanation. 

In order to understand the behavior of Sgs for low 
J we have calculated the profiles of the magnetization, 
(S{), in the spin sector S z = Sgs > and of the hole den¬ 
sity, ( rih,i ), for J = 0.1 (see Fig. [2]). To help in the data 
visualization, we use a linearized version of the lattice, 
as illustrated in Fig. Ha). As shown in Fig. H(b), for 
8 = 0.04 the holes distort the ferrimagnetic structure, 
which display a modulation with wavelength A ss 17, in 



FIG. 4. (Color online). Phase separation (IC-RVB) for a 
chain with N = 3 N c + 1 = 100 sites, J = 0.3, and Nh = 18 
holes: spin correlation function between B spins at the same 
cell, ( S Bl ,i ■ S b 2 ,i)i an d hole density profile, 


anti-phase with that exhibited by the hole (charge) den¬ 
sity wave. We have thus identified a modulated itinerant 
ferrimagnetic phase in this underdoped regime. On the 
other hand, as shown in Fig. He), for <5 = 0.18 the mag¬ 
netization has local maxima in coincidence with those of 
the hole density profile. In this case, the IC phase is char¬ 
acterized by the presence of ferromagnetic Nagaoka spin 
polarons [HI, [39 1 ] due to hole density wave with A « 4. 
Our results point to a value of J (~ 0.2) below which 
ferromagnetic “bubbles” appear as precursors of the F 
phase found for J < Jp,s [see Fig. Ha)]. 

For J = 0.3, Sgs = 0 in the IC phase, as shown in 
Fig. He). In Figs. [4] (a) and (b) we present the magnetic 
structure factor 


S(q) = 


1 


S L (S L + 1) 


2AT C + 1 

£ e iq( ' I ~ m \S l ■ S m ), (2) 


l,m 


where l, m and S refer to the lattice representation shown 
in Fig. H a )i f° r this value of J and doping ranging from 
8 = 0 up to 5 = 0.12. In a long-range ordered ferrimag¬ 
netic state, sharp maxima at q = 0 (ferromagnetism) and 
q = 7 r (antiferromagnetism) are observed in the curve 
S(q) for (5 = 0. Adding two holes to the undoped state, 
sharp maxima at q = 0 and tt are also observed, while 
broad maxima occur for 8 = 0.04, indicating short-range 
ferrimagnetic order which evolves to the IC phase by in¬ 
creasing doping, before phase separation (IC-RVB) at the 
line 8 = 8ps,j [see Fig. H a )]- In the inset of Fig. Hb) we 
show the departure of the maximum of S(q) from q = n. 


IV. PHASE SEPARATION, RVB STATES AND 
LUTTINGER LIQUID 

In Fig. Ha) the dashed line inside the PS region fix the 
boundary between two types of phase separation: in one 
case, the separation occurs between Nagaoka ferromag¬ 
netism and short-range RVB states (F-RVB); while in the 
other, it occurs between IC and short-range RVB states 
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FIG. 5. (Color online). Luttinger liquid behavior for a chain 
with N = 3N C = 24 (ED results), (a) Ratio R = u p /\JD x/f 
as a function of S for the indicated values of J. (b) Exponent 
K p as a function of 5. 


(IC-RVB). Indeed, for 0 < J < 0.063 and 8f-rvb < 
8 < 1/3, the GS phase separates with F and short range 
RVB states under coexistence, where Sf-rvb denotes 
hole density values along the phase separation line F- 
RVB, thereby extending our previous result fl3j valid 
only for J = 0. However, for 0.063 < J < 0.4 the sys¬ 
tem behaves differently. The new PS (IC-RVB) region is 
here illustrated for J = 0.3, N = 3 N c + 1 = 100 sites, 
and Nh = 18 holes: we thus find that there are 26 cells 
with odd parity (pi = —1), associated with the IC phase, 
and the remaining 7 cells with even parity (pi = +1), 
associated with the RVB phase. In this case, as shown 
in Fig. He), the hole-poor IC phase presents a local 
spin correlation function (S Pl ,i • Sb 2 ,i) ~ 0.2, average 
hole density per site ~ 0.16, estimated from the sites 
indicated by arrows [one A site and two B sites in the 
context of the effective linear chain shown in Fig. Ha)], 
and hole-density wave with A ~ 4; while the hole-rich 
RVB phase presents (S b±,i • S b 2 ,i) ~ —0.4 and average 
hole density per site ~ 1/3, estimated from a cell with A 
and B sites indicated by arrows. Therefore, apart from 
boundary effects, the above results thus indicate that the 
phase separation for a given J value is defined by the co¬ 
existence of the two phases with the hole densities Sic-ps 
(« 0.16 for J = 0.3) and S PS -rvb (« 1/3 for J = 0.3) 
fixed at the IC-PS and PS-RVB boundaries, respectively, 
while the size of the phases are fixed by the chemical 
doping 8 = Nh/N (= 0.18 for N = 100 and Nh = 18). 
We also remark that the stable RVB phase observed at 
<5 = 1/3 and 0 < J < 0.4, which has finite charge and 
spin gaps, is in agreement with predictions for J = 0.35 
PI and J = 0 [13. 

For 0 < J < 0.4 and 1/3 < 8 < 2/3, a crossover 
region with the presence of long-range RVB states after 
hole addition away from 8 = 1/3 is observed [see Fig. 
Ha)]- At the commensurate filling 5 = 2/3, the system 
presents a charge gap, while the spin excitation is gapless, 
also extending our previous result for J = 0 [l3| . 

With the aim of investigating the LL behavior as a 
function of J and 5 > 2/3, we have calculated, through 


0 . 03 ( b ) 



FIG. 6. (Color online), (a) Shift Egs~Ef per unit cell, where 
Ef is the energy of the fully polarized ferromagnetic state, as 
a function of J/Jf,s for the indicated values of N c and 6. (b) 
Instability line of the Nagaoka ferromagnetic phase, (c) Spin 
and (d) hole profiles, {rih : i) and (Si), respectively, for a chain 
with N = 3 N c + 1 = 100, 5 = 0.1, and the indicated values of 
S z and J. 


ED techniques, the ratio R = u p /yjD\/w, where 
N 

X = -f[E(N h + 2) + E(N h - 2) - 2 E GS (N h )} (3) 

is the charge susceptibility, and E(N} t ± 2) is the total 
energy for A/, ± 2 holes; 


N c \d 2 E(<S>) 
47T < 94> 2 


(4) 


is the Drude weight, where -E(4>) is the lowest energy for 
a magnetic flux $ through a closed chain, and <1 'min its 
value at Egs\ 


E(kcs + Afc, 5 = 0) — Ecsikcs, Sgs = 0) / ^ 

-TT- ( 5 ) 


is the charge excitation velocity, where A k = 2n/N c , and 
E(kcs + A k, S = 0) is the lowest energy with wavenum¬ 
ber k = kes + A k and total spin 5 = Sgs = 0. If 
the low-energy physics of the system is that of a LL, we 
should find R = 1 jjtQj; moreover, the exponent govern¬ 
ing the asymptotic behavior of the correlation functions, 
K p , satisfies the relation K p = nu p /2x ■ As shown in 
Fig. [^a) , R is indeed very close to 1 for 8 > 2/3; in 
addition, as shown in Fig. H^); we find 0.7 > K p > 0.5 
for 8 > 2/3. Remarkably, as shown in Figs. Ha) and (b), 
the data for R and K p exhibit data collapse as a function 
of 8 for 0 < J < 0.4. In short, the results above clearly 
indicate that for 8 > 2/3 and 0 < J < 0.4 the system 
behaves as a LL in the strong coupling regime. 
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V. STABILITY OF NAGAOKA 
FERROMAGNETISM 

In this Section, we shall provide strong evidence that 
for 0 < J < Jf,s and 0 < 8 < Sf-rvb-, the kinetic energy 
of holes is lowered in a fully polarized ferromagnetic state, 
an extension of Nagaoka ferromagnetism [23j, [24j], with 
the GS energy equal to that of non-interacting spinless 
fermions: Egs = Ef- 

The estimate of Jf,s is based on the data for the shift 
( Eqs ~ E f )/N c as a function of J, as illustrated in Fig. 
Ha) for 8 close or equal to 0.1. We stress that the shift 
decreases as N c increases for 0 < ( J/Jf,6 ) < 1, and goes 
to zero in the thermodynamic limit. In addition, one 
should notice that, by examining the data above and be¬ 
low J = Jf,s, particularly for N = 3 N c + 1 = 364 sites, 
dEcs/dJ appears to be discontinuous at J = Jf,s in the 
thermodynamic limit, thus suggesting a first-order tran¬ 
sition to the IC phase at ( J/Jf,s ) = 1- In Fig. Hb) we 
show that our estimated transition line, Jfs , [see also 
Fig. Ha)] is almost not affected by finite size effects and 
implies 8fj ~ yfj as 8 —> 0, as found from analytical 
results [23,123 for the t-J model in a square lattice. In 
particular, for J = 0 the instability of the Nagaoka state 
occurs at 5 w 0.23, which is very close to the values of 
hole doping estimated for n-leg ladder systems [25| and 
the square lattice snn. 

The spin profile for a chain with N = 3 N c + 1 = 100, 
J > 0 and 8 = 0.1 is also in very good agreement with the 
Nagaoka state, as shown in Fig. H c ); although boundary 


effects are visible for J > 0; in fact, S z changes from 45 
to 39 (on average, three spins at each boundary are not 
fully polarized), but one should notice that the change 
saturates as J slightly increases above zero. This fact 
is corroborated by the hole density shown in Fig. Hd), 
whose data for the referred states with S z = 39 at 8 = 0.1 
are very well described by the Nagaoka profile. 

VI. DISCUSSION AND CONCLUDING 
REMARKS 

The presented phase diagram of doped AB 2 t-J chains 
is remarkably rich. Indeed, several magnetic and non¬ 
magnetic phases manifest themselves in a succession of 
surprising relevant features, some of which are similar to 
those observed in the square and n-leg ladder lattices: 
all in a simple doped chain. In particular, we empha¬ 
size the modulated ferrimagnetic structure, the occur¬ 
rence of Nagaoka spin polarons in the underdoped regime 
and small values of J, incommensurate structures with 
nonzero magnetization, the strong-coupling LL physics 
in the high-doped regime, and the instability of Nagaoka 
ferromagnetism against J and doping. Therefore, these 
chains are unique systems and of relevance for the physics 
of polymeric materials, whose properties might also rep¬ 
resent challenging topics to be explored via analog simu¬ 
lations in ultracold fermionic optical lattices. 

This work was supported by CNPq and FACEPE 
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[1] The Hubbard Model - A Reprint Volume, edited by A. 
Montorsi (World Scientific, Singapore, 1992) 

[2] F. Gebhard, The Mott Metal-Insulator Transition 
(Springer, Berlin, 1997)M. Imada, A. Fujimori, and 
Y. Tokura, Rev. Mod. Phys. 70, 1039 (1998) 

[3] A. Auerbach, Interacting Electrons and Quantum Mag¬ 
netism (Springer, Berlin, 1998) 

[4] P. W. Anderson, Science 235, 1196 (1987); The Theory 
of Superconductivity in the High-Tc Cuprates (Princeton 
University Press, Princeton, 1997); P. A. Lee, N. Na- 
gaosa, and X.-G. Wen, Rev. Mod. Phys. 78, 17 (2006) 

[5] E. H. Lieb and F. Y. Wu, Phys. Rev. Lett. 20, 1445 
(1968)F. H. L. Essler, H. Frahm, F. Gohmann, A. 
Klumper, and V. E. Korepin, The One-Dimensional Hub¬ 
bard Model (Cambridge University Press, Cambridge, 
2005) 

[6] E. H. Lieb, in The Hubbard model, its Physics and Math¬ 
ematical Physics, Nato ASI, Series B: Physics, Vol. 343, 
edited by Baeriswyl, D. K. Campbell, J. M. P. Carmelo, 

F. Guinea, and E. Louis (Plenum, New York, 1995)H. 
Tasaki, J. Phys.: Condens. Matt. 10, 4353 (1998)G.-S. 
Tian, J. Stat. Phys. 116, 629 (2004) 

[7] E. H. Lieb, Phys. Rev. Lett. 62, 1201 (1989) 

[8] E. Lieb and D. Mattis, J. Math. Phys. 3, 749 (1962) 

[9] A. M. S. Macedo, M. C. dos Santos, M. D. Coutinho- 
Filho, and C. A. Macedo, Phys. Rev. Lett. 74, 1851 


(1995) 

[10] G.-S. Tian and T.-H. Lin, Phys. Rev. B 53, 8196 (1996) 

[11] S.-D. Liang, Z. D. Wang, Q. Wang, and S.-Q. Shen, 
Phys. Rev. B 59, 3321 (1999) 

[12] G. Sierra, M. A. Martin-Delgado, S. R. White, D. J. 
Scalapino, and J. Dukelsky, Phys. Rev. B 59, 7973 
(1999)M. A. Martin-Delgado, J. Rodriguez-Laguna, and 

G. Sierra, ibid. 72, 104435 (2005) 

[13] R. R. Montenegro-Filho and M. D. Coutinho-Filho, 
Phys. Rev. B 74, 125117 (2006)M. H. Oliveira, E. P. Ra- 
poso, and M. D. Coutinho-Filho, ibid. 80, 205119 (2009) 

[14] A. A. Lopes and R. G. Dias, Phys. Rev. B 84, 085124 
(2011)A. A. Lopes, B. A. Z. Antonio, and R. G. Dias, 
ibid. 89, 235418 (2014) 

[15] For Bose-Hubbard models, see J. J. Garca-Ripoll 
and J. K. Pachos, New Journal of Physics 9, 139 
(2007)S. Takayoshi, H. Katsura, N. Watanabe, and 

H. Aoki, Phys. Rev. A 88, 063613 (2013) 

[16] E. P. Raposo and M. D. Coutinho-Filho, Phys. Rev. Lett. 
78, 4853 (1997)Phys. Rev. B 59, 14384 (1999)C. Vitori- 
ano, M. D. Coutinho-Filho and E. P. Raposo, J. Phys. 
A: Math. Gen. 35, 9049 (2002) 

[17] F. C. Alcaraz and A. L. Malvezzi, J. Phys. A: Math. Gen. 
30, 767 (1997) 

[18] R. R. Montenegro-Filho and M. D. Coutinho-Filho, 
Physica A 357, 173 (2005)S. Yamamoto and J. Ohara, 



6 


Phys. Rev. B 76, 014409 (2007) 

[19] N. Ivanov, Condens. Matter Phys. 12 , 435 (2009) 

[20] R. R. Montenegro-Filho and M. D. Coutinho-Filho, 
Phys. Rev. B 78, 014418 (2008); K. Hida and K. 
Takano, ibid. 78, 064407 (2008); A. S. F. Tenorio, 
R. R. Montenegro-Filho, and M. D. Coutinho-Filho, 
ibid. 80, 054409 (2009)T. Shimokawa and H. Nakano, 
J. Phys. Soc. Japan 81, 084710 (2012)S. C. Furuya and 
T. Giamarchi, Phys. Rev. B 89, 205131 (2014) 

[21] M. S. S. Pereira, F. A. B. F. de Moura, and M. L. Lyra, 
Phys. Rev. B 77, 024402 (2008)79, 054427 (2009)0. Ro¬ 
jas, S. M. de Souza, and N. S. Ananikian, Phys. Rev. E 
85, 061123 (2012) 

[22] M. D. Coutinho-Filho, R. R. Montenegro-Filho, E. P. 
Raposo, C. Vitoriano, and M. H. Oliveira, J. Braz. Chem. 
Soc. 19, 232 (2008) 

[23] Y. Nagaoka, Phys. Rev. 147, 392 (1966) 

[24] H. Tasaki, Prog, of Theor. Phys. 99, 489 (1998) 

[25] L. Liu, H. Yao, E. Berg, S. R. White, and S. A. Kivelson, 
Phys. Rev. Lett. 108, 126406 (2012) 

[26] F. Becca and S. Sorella, Phys. Rev. Lett. 86, 3396 (2001) 

[27] E. Eisenberg, R. Berkovits, David A. Huse, and B. L. 
Altshuler, Phys. Rev. B 65, 134437 (2002) 

[28] M. M. Maska, M. Mierzejewski, E. A. Kochetov, L. Vid- 
mar, .J. Bonca, and O. P. Sushkov, Phys. Rev. B 85, 


245113 (2012) 

[29] E. Dagotto, Science 309, 257 (2005) 

[30] S. Chakrabarty, V. Dobrosavljevic, A. Seidel, and 
Z. Nussinov, Phys. Rev. E 86, 041132 (2012) 

[31] J. M. Tranquada, Physica B 407, 1771 (2012) 

[32] V. J. Emery, S. A. Kivelson, and H. Q. Lin, Phys. Rev. 
Lett. 64, 475 (1990) 

[33] P. Corboz, S. R. White, G. Vidal, and M. Troyer, 
Phys. Rev. B 84, 041108 (2011) 

[34] W.-J. Hu, F. Becca, and S. Sorella, Phys. Rev. B 85. 
081110 (2012) 

[35] P. Corboz, T. M. Rice, and M. Troyer, arXiv: 1402.2859 

[36] M. Ogata, M. U. Luchini, S. Sorella, and F. F. Assaad, 
Phys. Rev. Lett. 66, 2388 (1991) 

[37] S. R. White, Phys. Rev. B 48, 10345 (1993); U. Scholl- 
wock, Rev. Mod. Phys. 77, 259 (2005) 

[38] F. D. M. Haldane, J. Phys. C 14, 2585 (1981); J. Voit, 
Rep. Prog. Phys. 58, 977 (1995); T. Giamarchi, Quan¬ 
tum Physics in One Dimension (Oxford University Press, 
New York, 2003) 

[39] A. Auerbach and B. E. Larson, Phys. Rev. Lett. 66, 2262 
(1991); E. Dagotto and J. R. Schrieffer, Phys. Rev. B 43, 
8705 (1991) 

[40] See, e. g., C. A. Hayward and D. Poilblanc, Phys. Rev B 
53, 11721 (1996) 


